Tension between SNela and BAO: current status and future forecasts 
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Using real and synthetic Type la SNe (SNela) and baryon acoustic oscillations (BAO) data 
representing current observations forecasts, this paper investigates the tension between those probes 
in the dark energy equation of state (EoS) reconstruction considering the well known GPL model 
and Wang's low correlation reformulation. In particular, here we present simulations of BAO data 
from both the the radial and transverse directions. We also explore the influence of priors on ilm 
and Qb on the tension issue, by considering la deviations in either one or both of them. Our 
results indicate that for some priors there is no tension between a single dataset (either SNela 
or BAO) and their combination (SNela+BAO). Our criterion to discern the existence of tension 
(cr-distance) is also useful to establish which is the dataset with most constraining power; in this 
respect SNela and BAO data switch roles when current and future data are considered, as forecasts 
predict and spectacular quality improvement on BAO data. We also find that the results on the 
tension are blind to the way the GPL model is addressed: there is a perfect match between the 
original formulation and that by the correlation optimized proposed in [27], but the errors on the 
parameters are much narrower in all cases of our exhaustive exploration, thus serving the purpose 
of stressing the convenience of this reparametrization. 
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I. INTRODUCTION 

Almost thirteen years ago, the accelerated expansion 
of the Universe was discovered while reconstructing the 
Hubble diagram of SNeIa[T]. From then on, a large 
amount of cosmological data has been collected: the 
already mentioned Hubble diagram of SNela the 
measurements of cluster properties as the mass, the 
correlation function and the evolution with redshift of 
their abundance |3]; the optical surveys of large scale 
structure [4]; the anisotropics in the cosmic microwave 
background (CMB) |3 E]; the cosmic shear measured 
from weak lensing surveys [7] and the Lyman - a forest 
absorption [5]. All these data sets confirm that the 
Universe is spatially flat, exhibits a subcritical matter 
content and its expansion is accelerated; but despite all 
these recent progresses, the origin of this accelerated 
expansion remains unknown. 

The simplest and most accepted model is the ACDM 
one, where the acceleration is driven by a positive cosmo- 
logical constant, A [SKTU], which has to be small enough 
for it to start dominating the Universe only at late times, 
where it represents about the 73% of the total energy con- 
tent of the Universe, as reported by the WMAP7-year 
analysis [TT]. At the same time, it requires the presence 
of a large amount of cold dark matter (about the 23% 
of the total matter-energy content), non-baryonic matter 



which is detectable only by its gravitational interaction 
with ordinary baryonic matter. The ACDM model pro- 
vides a good fit to most of the data [H] , and still remains 
the best candidate in some respects, as demonstrated by 
the WMAP7-year data [B]. But it is also affected by seri- 
ous theoretical shortcomings [131 [H] that have motivated 
the search for alternative candidates generically referred 
to as dark energy. Attempts in this direction include, 
among others, modified gravity |15| . Cardassian cosmol- 
ogy [T^ , the Chaplying gas [T7] , braneworld models [TH] , 
and f{R) theories [15] . 

Studying the nature and the main properties of dark 
energy practically means studying its equation of state 
(EoS), w{a) = px{a)/pxia), where px and px are re- 
spectively its density and pressure and a = 1/(1-1- z) is 
the scale factor. It follows then that the dark energy 
density function (in critical density units) is 



fix (a) oc exp 



da' 



which appears in the Hubble function as 

H{a) = Ho [flraa-^ + nx{a)\ 



1/2 



(1) 



(2) 
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where f2„i is the present value of the matter density 
in critical density units and (1 — ilm) = ^x{o- = 1) 
is the current dark energy density in the same units. 
Clearly, the spatial flatness hypothesis has been made, 
and radiation and curvature contributions have been 
neglected. 

This last equation makes it possible to study and con- 
strain the dark energy EoS within a suitable model and 
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relate it to observations, although actually extra ingredi- 
ents such as the fraction of baryons fif, or the comoving 
sound horizon size may be needed. The most important 
and difficult goal at this stage is to address a physically 
adequate expression for the dark energy EoS. To be use- 
ful, w(z) must be sufRciently sophisticated to be able 
to explain the data, and simple enough so as to provide 
reliable predictions. Eventually, once the free parame- 
ters in the chosen dark energy model and the remaining 
cosmological parameters have been observationally con- 
strained, a picture of the evolution of the Universe will 
emerge. 

It is quite strongly stablished that dark energy domina- 
tion began somewhat recently, and therefore low redshift 
data, are precisely those best suited for its analysis. The 
two main astrophysical tools of such nature are standard 
candles (objects with well determined intrinsic luminos- 
ity) and standard rulers (objects with well determined 
comoving size). Such probes provide us with distance 
measures related to H{z), and the best so far represen- 
tatives of those two classes are SNela and BAO. Those 
are in fact low redshift datasets, and much effort is being 
done in those two observational contexts toward obtain- 
ing more and better measurements [20] . 

It is well known, however, that a "tension" pTl 
among the estimated EoS parameters from different 
datasets in general and those two mentioned in partic- 
ular can arise. The word "tension" will be used in this 
paper in its more general meaning, namely, to indicate 
that the EoS parameters values obtained by using a cer- 
tain dataset can differ from those obtained from another 
dataset at least at a 2a level. Typically, BAO data seem 
to prefer phantom dark energy [w < —1), which in fact 
is the key point to be addressed, as a confirmation would 
definitely do away with the ACDM model. On the other 
hand, SNela data seem to favour a larger proportion of 
dark matter in the Universe than BAO data (which un- 
like supernovae are both connected to the cosmological 
background and also to its structure). Clearly, tension 
between those datasets becomes manifest in rather rele- 
vant aspects of dark energy reconstruction. 

As long as tension is concerned, one cannot waive some 
additional facts. First of all, priors in the fraction of dark 
matter and/or baryons may influence the estimated val- 
ues of the dark energy parameters rather strongly , as dis- 
cussed in [ini US] . Secondly, the tension may be bringing 
to light some caveats in the dark energy parametrization. 
Thirdly, it may just be a statistical pathology produced 
by the limitations of the observational data available [23] . 

All those questions arise ,and hopefully will be given 
some answers, in the context of any of the particular 
parametric forms of •w{z) [24j (just to cite a few) which 
have been proposed to describe the behavior of time- 
evolving dark energy EoS. At present, due to the cur- 
rent limited constraining power of the data not one of 
those parametrizations has championed definitively and 
there is still a long of room for discussion about their 
suitability either along the focus or many other inter- 



esting routes. Nevertheless, practicality has made some 
models particularly popular |25| , with the the Chevallier- 
Polarski-Linder (CPL) model [T31 US] excelling among 
them. It is given by 

w{a) = wq + {1- a)wa , (3) 

where wq is the present value of the EoS and Wa con- 
straints the value of the EoS at early times, i.e, for a — 
we have w ^ Wg + Wa- This scenario has gathered many 
followers because it is simple and sensitive enough to 
good quality astronomical data. This parametrization 
is consistent with low redshift data as SNela and BAO 
data, although some unusual effects manifest themselves 
when CMB data are considered (in the vicinity of the 
regions of the parameter space where dark energy mim- 
icks dark matter). Nevertheless, Wq and Wa are highly 
correlated, thus making it harder to extract very strong 
conclusions about dark energy evolution from them. 

A reformulation of the CPL model was proposed in 
[27] which puts the accent on the evaluation of the same 
dark energy EoS at two different redshifts, one being the 
present one wq = w{a = 1) = w{z = 0) and the other 
being that corresponding to a pivot redshift value: Wc = 
w{z = Zc)- This pivot value is chosen appropriately so 
that the correlation between wq and Wc is minimal, as 
this turns to be advantageous. 

Motivated by the issue of the tension among SNela and 
BAO data, and taking into account the points already 
made above in this section, we work along the following 
main directions. First, we study how the choice of pri- 
ors affects the basic features of our dark energy models 
(separation from the ACDM model). Then, we explore 
whether the influence of the choice of priors gets atten- 
uated when the correlation between parameters is mino- 
rated. After that, we repeat those tests with synthetic 
data based on the speciflcations of future survey proper- 
ties, thus checking the influence of future systematics on 
EoS estimations and errors. The mock data we present 
are to the best of our knowledge the flrst mock collections 
of data reproducing future observations derived from the 
radial and transverse BAO modes. 

The paper is organized as follows; in Sec.|ll]we briefly 
review the EoS parametrizations considered; in Sec. |III| - 
IIVI we describe the observational and simulated mock 
datasets; in S ec. [V] we review the used statistical method- 
ology; in Sec. |VI| we outline and comment the results of 
our analysis. 

II. LOW CORRELATION DARK ENERGY EOS 

In [27] the popular CPL parametrization was given an- 
other turn of the screw to become 

= (^) + (^) 

where Wq = w{z — 0) and Wc = w{z = Zc). The moti- 
vation for this reformulation, as we have already pointed 
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out and will develop further later, is to decorrelate the 
parameters as much as possible. Taking into account 
Eq.Q, and assuming the existence of a matter compo- 
nent as well, the square of E{z) = H{z)/Ho reads 



E^z) = n„,{i + zf + (1 - n„,){i + z)3[i+(™T^ 

, o l^'^c -■Wo\ Z 

X exp 3 



flc — 1 J 1 + Z 



(5) 



Eq.Q can be related to the usual CPL parametrization 

by 



(6) 



In principle, the subindex c should indicate the scale fac- 
tor (or redshift) value for which the parameters (wq, Wc) 
are uncorrelated. However, this value depends on the 
single dataset (or combination of them) one is consider- 
ing; in 1^7] the author proposed to fix it as z^ — 0.5, 
this value being sufficiently close to the one derived from 
current data {zc ^ 0.3) and thus arguing that the corre- 
lation between {wqtWc) will consequentially be relatively 
small. Moreover, starting from the definition of Wc, it is 
straightforward to demonstrate that the couple (wq^Wc) 
will always be less correlated than the CPL parameter 
couple {wQ,Wa), when 



(T^(wo) < 2|(1 - ac)a'^{woWa)\- 



(7) 



Note that this condition is always satisfied for Zc = 0.5, 
i.e. flc = 2/3. For this particular value of Oc one gets 



E^{z) = n„r{l + Z)^ + (1 - a„)(l + 2)3(l-2t«o+3«.o.5) 
9{wo - Wo.5)z 



X exp 



1 + z 



(8) 



Before closing let us recall that we will at all stages be 
working with two testbenches: Wang's (low correlation) 
model on the one hand, and the CPL model on the other 
(we omit the specifics of E{z) for the CPL model as this 
has been reproduced in the literature extensively). 



III. CURRENT ASTROPHYSICAL DATA 

For what concerns the kind of observational data we 
are going to use, we remind that the most familiar and 
sensitive probes for testing dark energy properties are the 
SNela: even though they are extremely rare astrophysi- 
cal events, the modern and specifically planned strategies 
of detection make it possible to observe and collect them 
up to relatively high redshifts (z ~ 2) well beyond 
the era when dark energy should have become the dom- 
inant energy term in the Universe. Another important 
property of modern surveys is that they offer a statis- 
tically significant number of observations, whose quality 
has improved considerably over the years [29) . Present 
samples are made of approximately 500 events; future 
surveys are planned to detect up to approximately 5000 



events. However, these SNela datasets are not always 
consistent with other types of cosmological observations 
(as we shall see), and even an historical tension can be 
detected among different SNela samples [HI [55] . 

While SNela are well known to be standard candles, 
we also have standard rulers; among these, one of the 
main techniques rests on the BAO peaks detection in 
the galaxy power spectrum. This has recently emerged 
as a promising standard ruler for cosmology, potentially 
enabling precise measurements of the dark energy pa- 
rameters with a minimum of systematic errors [30) 131) . 
For this technique to be useful, one must explore, as uni- 
formly as possible, a large volume in real space with a 
sufficiently high sampling density. These requirements 
are met by current high quality data (the Sloan Digi- 
tal Sky Survey (SDSS) [51 [33] and the 2-degree Field 
Galaxy Redshift Survey (2dFGRS) (M]), and thus BAO 
have become a nature tool to study the cosmic expansion 
history. 



A. Supernovae 

We use the most recent Snela sample available, the 
Union2 sample described in [3S]. The Union2 SNela 
compilation is the result of a new dataset of low-redshift 
nearby-Hubble-flow Type la SNe and is built with new 
analysis procedures to work with several heterogeneous 
SNela compilations. It includes the Union data set from 
[36] with six added SNela first presented in [35], along 
with SNela from [37], the low-z and the intermediate-z 
data from [3H] and [32] respectively. After the application 
of various selection cuts to create a homogeneous and 
high signal-to-noise data set, A/sn = 557 SNela events 
distributed over the redshift interval 0.015 < z < 1.4 
were obtained. 

The statistical analysis of the Union2 SNela sample 
rests on the definition of the modulus distance: 



m(2j,Mo) = 51ogioML(2j, 1^™;^)] +Mo, 



(9) 



where dL{zj,flm,',d) is the Hubble free luminosity dis- 
tance: 

d.(.,a„;e)Hi + ^)f d.'^^^-^. (10) 

With the chosen notation we clarify the different roles 
of the various cosmological parameters appearing in the 
formulae: the matter density parameter appears sep- 
arated as it is assumed to be fixed to a prior value, while 6 
is the EoS parameters vector {0 — (1^0,^0.5) for Wang's 
model and = {wq, Wa) for the CPL one), which are the 
parameters that we will be constraining. The best fits 
will be obtained by minimizing the quantity 



Xsn(mo,6') 



L 



{ll(Zj,Vtm\ Ho, 9)} - y.obs{Z])f 



(11) 
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where the af, , are the measurement variances. The nui- 
sance parameter /xq encodes the Hubble parameter and 
the absolute magnitude M, and has to be marginalized 
over. Giving the heterogeneous origin of the Union2 
dataset, and the procedures to reduce data described 
in |36|, we will work with an alternative version [40] of 
Eq. (11), which consists in minimizing the quantity 



C3 



with respect to the other parameters. Here 



Cl 



C2 = 



E 
E 



{n{Zj,flm]fJ-0 = 0,0)} - ^lobs{Zj)Y 



{fl{Zj,n,n;fJ.O = 0,0)} - Hobs{Zj)) 



C3 



E ^ 



(12) 



(13) 



(14) 



(15) 



It is easy to see that Xsn just a version of Xsnj min- 
imized with respect to iiq. To that end, it suffices to 
notice that 



Xsn(Mo,0) = Cl - 2c2^o 



C3Mo 



(16) 



which clearly becomes minimum for /ip = C2/C3, and so 
we can see Xsn = Xsn(A''0 — 0,9). Furthermore, one 
can check that the difference between Xsn ^^'^ Xsn 
negligible. 



B. Baryon Acoustic Oscillations 

BAO have rarely been used independently of other 
datasets (such as SNela or the CMB) due to the small 
number of datapoints available so far |41j . but they pro- 
vide a measuring stick for a better understanding of 
the nature of the accelerated expansion of the Universe. 
Moreover, they stand on a completely independent basis, 
as compared to the SNela technique, and can contribute 
important features by comparing the data of the sound 
horizon today (using the galaxies clustering) to the sound 
horizon at the time of recombination (extracted from the 
CMB). 

In [42] the authors analyze the clustering of galax- 
ies within the spectroscopic Sloan Digital Sky Survey 
(SDSS) Data Release 7 (DR7) galaxy sample, including 
both the Luminous Red Galaxy (LRG) and Main sam- 
ples, and also the 2-degree Field Galaxy Redshift Survey 
(2dFGRS) data. In total, the sample comprises 893319 
galaxies over 9100 deg2. For a redshift survey in a thin 



shell, the position of the BAO peak approximately con- 
strains the ratio 



rsjzd) 
Dv{z) 



(17) 



where rs{zd) is the comoving sound horizon at the baryon 
dragging epoch 



Azd) 



Ho 



Csjz) 

E{z) 



dz 



(18) 



with c the light velocity, Cg the sound speed and Zd the 
dragging epoch redshift. By definition, the dilation scale 
Dv{z) is 



Dv{z,Vl,n;9) = 



(1 



z?Dl 



^H{z,^m]e) 
where Da is the angular diameter distance 



1/3 



DA{z,n„,-0) = 



1 



l + z 



H{z',vi^-e) 



(19) 



(20) 



Again in this case we will choose fi™ to be fixed to the 
chosen prior value, while 6 will be the EoS parameters 
vector. Through the comoving sound horizon, the dis- 
tance ratio dz is related to the expansion parameter h 
(defined such that H = lOO/i) and the physical densities 
ilm and ilf,. Specifically, following [351135] we have 



rs{zd) = 153.5 



0.02273 



-0.134 



0.1326 ) 



-0.255 



Mpc 



(21) 

Modeling the comoving distance-redshift relation as a cu- 
bic spline in the parameter Dv{z), in [32] the authors 
obtained the following best fit results: 



do.2 = 0.1905 ± 0.0061, 
do.35 = 0.1097 ± 0.0036, 

with an inverse covariance matrix 



^BAO 



30124 -17227 
-17227 86977 



The function for BAO is defined as 

Xbao(^) = ^bao(^)Cbao^bao(^), 
where Xbao is defined as 



X 



BAO 



'•s(zd) 

i3v-(o.2,n„;e) 

rs(zd) 
Dv(0.35,n„;e) 



0.1905 
' 0.1097 



(22) 
(23) 

(24) 
(25) 



IV. MOCK DATA 

Even though some astrophysical data like supernovae 
have a considerable quality and samples are made of a 
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significant number of datapoints, the derived constraints 
on dark energy are still not narrow enough, and therefore 
uncertainties on its nature and features are still consid- 
erable. This is particularly important as regards BAO 
observations, which is becoming a very active area with 
many planned future surveys [20]. Thus, forecasting the 
improvements to be expected from future missions is of 
paramount interest. The argued relevance of the ques- 
tions we propose in this paper suggests the usefulness 
of reexamining them in the light of synthetic data gen- 
erated under the specifications of future surverys. Our 
main tool for this task is the free software for cosmolog- 
ical calculations Initiative for Cosmology (iCosmo) (see 
[? ]) and its BAO modules. This package is the ofspring 
of considerable theoretical efforts in different forecast as- 
pects |31j , and we have modified and extended it to suit 
it to our needs. 



A. SNela simulations 

To create SNela mock samples we have taken into ac- 
count the specifications of the Wide-Field Infrared Sur- 
vey Telescope (WFIRStQ a space mission planned by 
NASA which has among its primary objective is to ex- 
plore the nature of dark energy employing three distinct 
techniques, measurements of BAO, SNela distances, and 
weak gravitational lensing. In [47] and Table |l] we report 
the redshift distribution of the expected observed SNela 
within various redshift bins. The main redshift range ex- 
tends up to z ^ 1.6, with 4 expected SNela at z > 1.6. 
The very low redshift subsample, i.e. the 500 SNela at 
z < 0.1, is supposed to come from the Nearby Super- 
nova Factory (SNfactory) [35] with measurements in the 
redshift range [0.03,0.08] so we are going to assume the 
former as the lowest redshift value in our mock sample. 

Our formula for errors on SNela magnitudes stems 
from a recipe used in the binning approach |49j . which 
we have adapted to the case with one supernovae per 
redshift bin: 



syst ' 



(26) 



where 



• CTrnt — 0.15 is the intrinsic dispersion in magnitude 
per SN, assumed to be constant and independent 
of redshift for all well-designed surveys; 

• Cpec = 5cr„/ (In(lO)cz) is the error due to the uncer- 
tainty in the SN peculiar velocity, with = 500 
km/s, c is the velocity of light and z is the redshift 
for any SN; 

• f^syst = 0.02(z/zmaa:) is the floor Uncertainty re- 
lated to all the irreducible systematic errors with 



cannot be reduced statistically by increasing the 
number of observations. The value 0.02 is con- 
servative from the perspective of what space-based 
missions could achieve. Those are precisely the re- 
sources expected to provide high-redshift SN, which 
are in turn the ones in which the systematic er- 
ror term is expected to contribute. Note as well 
that Zmax is the maximum observable redshift in 
the considered mission and this linear term in red- 
shift is used to account for the dependence with 
redshift of many of the possible systematic error 
sources (for example the Malmquist bias or gravi- 
tational lensing effects). 

In addition, we have included some extra (though very 
slight) noise, and then we have checked that our mock 
data are compliant with the main features and trends of 
the Union2 sample. 

For our analysis we have use three different fiducial 
models, all derived from the WMAP7-year analysis |50j : 
one is the quiessence model derived by using only CMB 
data and is named CMB-oriented model, another one is 
the quiessence model coming out from combining CMB 
data with BAO and is named BAO-oriented model, and 
the last one is the quiessence model coming out from 
combining CMB data with SNela and is named SNela- 
oriented model. The main difference among them is the 
value of the EoS parameter, w: in the CMB-oriented 
case it clearly corresponds to a phantom model, i.e. 
w < —1; in the BAO-oriented model the regime is phan- 
tom (though only slightly) , whereas in the SNela-oriented 
model w > —1. 

The function for the SNela mock data has been 



constructed as described in Sec. (Ill A). The supernovae 



^ http://wfirst.gsfc.nasa.gov/ 



redshift distribution we have follow is tabulated in (fol- 
lowing [5D]); then, through a simple extension to iCosmo 
of our own we managed to generate the values for 
those redshifts as drawn from hundreds of universes ran- 
domly and normally distributed around the fiducial one, 
and finally a reduction was been performed to give our 
mock sample. We again produced three mock samples 
using the previously mentioned fiducial scenarios. 

The use of different datasets offers a wider picture of 
the problem and allows for more robustness in our con- 
clusions, which on the other are only indicative and not 
as substantial as those coming from real data. 



B. BAO simulations 

Future BAO surveys will target at larger volumes and 
improved statistical features, so the result will be nar- 
rower constraints on dark energy. The current approach 
(as described in Sec. IIIB[ ) considers the quantity £>y, 
which is (modulo some constants) the geometric mean of 
the following two quantities which are becoming the new 
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Table I. Redshift distribution of WFIRST SNela 
samples 



Redshift bin # of SNela 



< 0.1 


500 


0.1 - 0.2 


200 


0.2-0.3 


320 


0.3-0.4 


445 


0.4-0.5 


580 


0.5-0.6 


660 


0.6-0.7 


700 


0.7-0.8 


670 


0.8-0.9 


110 


0.9- 1.0 


80 


1.0 - 1.1 


25 


1.1 - 1.2 


16 


1.2 - 1.3 


16 


1.3 - 1.4 


4 


1.4- 1.5 


4 


1.5- 1.6 


4 


> 1.6 


4 



observational targets 



rs{Zr) 



and 



rs{Zr) ' 



(27) 



where 



• r(z) is the comoving distance to redshift slice z and 
it is related to the transverse BAO mode; 

• r'{z) = dr(z)/dz = c/ {HqE{z)) is the derivative of 
r(z), and it is related to the radial BAO mode; and 

• rs{zr) is the sound horizon at recombination. 

In general, BAO measurements are more sensitive to dark 
energy evolution as they involve 1/H{z) directly, whereas 
luminosity data from SNela are related to an integral of 
that quantity so sensitivity will typically be somewhat 
compromised. 

The publicly available code has built-in routines based 
on the universal BAO fitting formulae for the diagonal 
errors on y and y' presented in [5T] . We have made the 
proper modifications of the code to replace the defaults 
with the EUCLID [51] survey properties as specified in: 
this survey is expected to cover 20000 sq. deg. and ob- 
serve 6.1 X 10^ galaxies in the redshift range 0.5 < z < 2.1, 
and as it is a spectroscopic survey we fix the redshift pre- 
cision as 0.001(1 -I- z). Then, we have written extra codes 
to generate a large number of normal random realizations 
around the mentioned three fiducial models Finally, af- 
ter some reduction, the syn t heti c BAO dat_a presented 
in Tables in Tables |lT]- ~ 
obtained. 



Ill - IV and Fig. ([T]) have been 



Table II. Simulated BAO data for y and y': CMB-oriented sam- 
ple. 



y 



y 



0.6 
0.8 
1.0 
1.2 
1.4 
1.6 
1.8 
2.0 



15.005 
19.221 
22.708 
25.762 
28.635 
31.126 
33.540 
35.357 



0.117 
0.109 
0.102 
0.098 
0.096 
0.101 
0.107 
0.112 



21.401 
18.953 
17.012 
15.098 
13.421 
12.213 
10.987 
10.094 



0.291 
0.186 
0.137 
0.101 
0.079 
0.069 
0.061 
0.055 



Table III. Simulated BAO data for y and y': BAO-oriented 
sample. 



z 

0.6 
0.8 
1.0 
1.2 
1.4 
1.6 
1.8 
2.0 



y 

15.021 
19.094 
22.768 
25.910 
28.751 
31.354 
33.725 
35.840 



0.118 
0.108 
0.103 
0.099 
0.096 
0.102 
0.107 
0.112 



y 

21.488 
19.147 
17.049 
15.216 
13.672 
12.340 
11.130 
10.154 



0.294 
0.189 
0.134 
0.101 
0.080 
0.070 
0.061 
0.055 



To define the function for BAO mock data we have 
to take into account that y and y' are correlated: such 
correlation can be quantified by a correlation coefficient 
p = (^yyijoyOyi « 0.4 [33]. With this in mind, the 
function is now defined as |44] 



Xbao(^) = 



1 



\2 



{y'{zj,n^,n^-e)~y'{z^)Y 



y'{zj,nm,^b;d) ~ y'{zj) 



(28) 



which clearly reduces to the more familiar expression for 
when the observable quantities are uncorrelated (i.e. 
P = 0). 



V. TWO-DIMENSIONAL GRID ANALYSIS 

Given the appropriate x^(^) function we then proceed 
to find the values of the model parameters which min- 
imize it. Using a grid (following prescriptions in |53j ) 
we study the parameter space of the model searching for 
the minimum. It is well known that when the num- 
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Table IV. Simulated BAO data for y and y': SNela-oriented Table V. Planck priors and CPL - Wang parametrization. Min- 
sample. imum correlation redshift Zc- 



z 


y 


ay 


y 


ay, 


0.6 


14.859 


0.116 


20.881 


0.284 


0.8 


18.770 


0.106 


18.624 


0.184 


1.0 


22.300 


0.101 


16.619 


0.131 


1.2 


25.442 


0.097 


14.854 


0.099 


1.4 


28.296 


0.095 


13.364 


0.078 


1.6 


30.810 


0.100 


12.089 


0.068 


1.8 


33.131 


0.105 


10.939 


0.060 


2.0 


35.181 


0.110 


9.983 


0.054 



0.0070 



0.0050- 



0.0030 




0.6 0.8 1.0 



1.2 1.4 
z 



1.6 1.8 2.0 



0.015- 




0.005- 



0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



Figure 1. BAO mock data: fractional errors of BAO quanti- 
ties y{z) and y'{z) for the Euclid survey. 



ber of parameters gets large, grid techniques require an 
increasing computational and time consuming effort, so 
it becomes preferable to use other statistical methods, 
such as Markov Chain Monte Carlo (MCMC). As we are 
going to explore two dark energy parameters only (with 
ilm and fib fixed), the grid method is suitable for our 
purposes. 

Combining some intuitive guesses with the knowledge 
provided by the vaste literature on the topic, we typi- 
cally require just two grids (the second being finer and 
narrower) to find a suitable initial point to feed the 
Levenberg-Marquadt algorithm (described in full detail 
in [S3], pp. 683 — 685), and eventually detect the min- 
imum value. In addition, the algorithm calculates 



Priors 




Zc (CPL/Wang) 


(f2m, ^b) 


SNela 


BAO 


SNela-l-BAO 


(0.214,0.0444) 


0.226 


0.113 


0.220 


(0.237,0.0444) 


0.222 


0.112 


0.215 


(0.260,0.0444) 


0.217 


0.112 


0.210 


(0.214,0.0405) 


0.227 


0.113 


0.220 


(0.237,0.0405) 


0.222 


0.112 


0.216 


(0.260,0.0405) 


0.217 


0.111 


0.209 


(0.214,0.0366) 


0.226 


0.112 


0.221 


(0.237,0.0366) 


0.222 


0.112 


0.216 


(0.260,0.0366) 


0.217 


0.111 


0.210 



the estimated covariance matrix C for the model param- 
eters 0, which we use to plot the confidence regions in 
the 2-dimensional parameter space under the Gaussian 
approximation. This means that provided 

Ax'(0) = {e- 0best)^C-l(0 - ©best), (29) 

the la (68.3%) and 2a (95.4% confidence levels) will cor- 
respond to Ax^ = 2.30 and Ax^ = 6.17 respectively. 

Finally, to compare results and test the issues of in- 
terest (in particular the tension among datasets), we 
compute the so called cr-distance, d^, i.e. the distance 
in units of a between the best fit points of the total 
sample and the best fit points of (respectively) SNela, 
BAO and a ACDM model corresponding to the different 
(rimjf^b) pairs that define our priors. Following [53^, the 
(T-distance is calculated by solving 

1 - r(l, |Ax^/2|)/r(l) = erHdjV2). (30) 

For homogeneity and consistency our "ruler" is in ev- 
ery case the total function (x^t = Xsn + a|ao)' 
and our recipe is the following: if we want to calcu- 
late the tension between SNela and SNeIa-|-BAO and 
the best fit parameter vectors are respectively 0sn and 
^SN-i-BAO then the previous Ax^ will be defined as 
Xtot(^SN+BAo) - Xtoti^SN)- Other cases follow simply 
from this prescription. 



VI. RESULTS AND CONCLUSIONS 

As we have discussed in Sec. |lj there are arguments 
suggesting that considering dark energy constraints from 
the combination of SNela and BAO data only is both 
relevant and useful, and so is comparing the individual 
predictions drawn from each other. This does not mean 
that we are going to avoid completely possible clues from 
the CMB analysis; in particular, the chosen priors for the 
cosmological density parameters flm and fib are derived 
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from a forecast of CMB observations with the Planck 
satelUte mission as reported in [54]. The predicted best 
fits and la errors are n^h'^ = 0.1308±0.0008 and f^b/i^ = 
0.0223 which, combined with our choice for h — 0.742 ± 
0.036 from ^J, become fi„ = 0.237 ± 0.023 and = 
0.0405 ± 0.0039. We then consider those best fits and 
combine them with their la errors errors to produce nine 
{ilm,^b) couples. We then feed these priors into our 
pipeline so that we can describe quite in detail how the 
estimated EoS parameters depend on the uncertainties 
in Final results for current data are shown 

in Tables ^VJ^ - Nuh, respectively for the CPL and 



Wang models; whereas forecasts for synthetic data can 
be examined by looking at Tables. (Villi - (IX I. 



Let us remind that Wang's EoS model is based on the 
idea of obtaining a minimum correlation among its pa- 
rameters. The pivot value wq.s is a conservative choice 
put forward in |27j . which achieves quite a low degree 
of correlation and provides simple expression, although 
each single dataset or combination considered will have 
its own optimal correlation redshift which typically will 
not be too far away from Wang's choice. The specific 
value of this pivot redshift can be easily calculated ap- 
plying error propagation to the definition of Wc in Eq. ^ , 
and then computing the correlation parameters between 
Wq and Wc as a function of Uc (or z^) which one finally 
minimizes. For the CPL model the final expression is |27] 



1 



a'^iwoWa) ' 



(31) 



which clearly depends on the data set considered. For 
real data (see Tab. ( VI l), we effectively we detect a mean 
value of Zc « 0.22 — 0.23 using our chosen SNela data set, 
which gets only slightly reduced to Zc ~ 0.21 — 0.22 when 
considering the total sample SNeIa-|-BAO, and ~ 0.11 
for BAO data. 



A. Tension: results 

The tension is confirmed both in the CPL and in the 
Wang case, and both with current and mock data; but 
it also shows different features with respect to references 



1. Current data 

In these cases, SNela and BAO best fits can differ by 
more than 2cr, as it is possible to conclude just by vi- 
sual inspection from Figs. ([2]) - (|3|. When we compare 
the best fits from SNela alone and the SNeIa-(-BAO com- 
bination case we conclude there is no tension, as in all 
cases we have a separation smaller than 2a. On the con- 
trary, there is always tension between BAO alone and 
SNeIa-|-BAO. Resorting again to visual examination we 
infer that for current data, the SNeIa-|-BAO credible in- 
tervals are very much infiuenced both in position and size 



by the SNela data, as these have much more statistical 
weight. 

We can also underline another important property: the 
tension depends strongly on the value of the cosmological 
priors. In particular, we find a precise trend: the a- 
distance between SNela and SNeIa-|-BAO best fits gets 
reduced either when or fif, grow. In other words, 
compliance between SNela and BAO data improves when 
the amount of dark energy decreases, as has been noticed 
in earlier works If and fif, are large enough 

(17™ = 0.260 and Qb = 0.0444, 0.0405), we even get close 
encounters (i.e. the distance among the best fits is less 
than la). The same trend holds true between SNela and 
BAO in the sense that the larger either ri„i or ilf, become, 
the larger the overlap between credible intervals. 

As expected, there is large sensitivity to changes in 
(from which both SNela and BAO depend) as compared 
to variations in Qb (which only influences BAO aspects) . 
Note that BAO contours are much larger in size thus 
reflecting their statistically less meaningful character. 



2. Mock future data 

When looking at mock data in Figs. ^ - (9]), we get 
tension again, but now the role of SNela and BAO is in- 
verted: the tension between SNela and SNela+BAO is 
always larger than 2a while the tension between BAO 
and SNeIa-|-BAO is much less than la when has the 
largest and 17 {, has the two largest values. 
Moreover, in this case the joint is clearly driven by 
BAO, exactly the opposite of the real data case. A pos- 
sible explanation for this is that, as discussed in the re- 
lated section, our mock data have been generated using 
different fiducial models, and two of them (the CMB and 
BAO-oriented case) have a constant value of w, below 
— 1, even if with large errors). This choice clearly cre- 
ates a bias in agreement with previous analyses: the fact 
that BAO are much keener with a phantom EoS than 
SNela, combined with the dramatically reduced errors 
drawn from expectations produce much smaller credible 
contours from BAO, very close in position and size to the 
total ones. In comparison, the statistical role of SNela is 
expected to become subdominant for future data. 

On the other hand, we have to note that we have a 
BAO driven joint even in the SNela-oriented case, 
where w > —1. This may be put down to the fact that 
even in this case the smearing effect in the luminosity 
data is strong enough to make the data less informative 
than the BAO ones. The clear tension between SNela 
and BAO is somewhat surprising given that we have gen- 
erated both synthetic datasets using exactly the same 
fiducial models. The element to have into consideration 
here is that the BAO data are divided in transverse and 
radial modes. The transverse modes share with SNela 
one fundamental property: they inform us on a quantity 
which depends on w{z) through an integral (as can be 
seen in Eq. (33)), and this makes them less sensitive to 
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changes in the EoS through the so-called smearing effect 
produced by the integration. Moreover the realistic dis- 
persion of data around the fiducial model induces some 
dissociation from the underlying model. On the other 
hand, the radial modes depend directly on w{z), so the 
BAO radial component can provide much more refined 
information. 

Assuming the possibility of a substantial smearing ef- 
fect, and considering the different statistical weights of 
SNela and BAO, one can wonder if there could possibly 
be a physical motivation for the tension among BAO and 
SNela, perhaps due to something not well established in 
the BAO or in the SNela data (for the SNela case, see 
[57] for a possible solution). 

Finally, looking at the cr— distances we can also verify 
that the SNela-oriented model seems to be somewhat dis- 
favored with respect to more phantom models, showing 
a larger tension (larger values for the cr— distances) with 
respect the other two cases. 



B. Dark energy evolution: results 

There is a number of interesting conclusions about 
the EoS parameters which one can draw from a general 
overview of our figures and tables. 

1. Pivot choice and effect on errors 

First of all, as Wang's model is a reparametrization of 
the CPL model, it was to be expected that many similari- 
ties between them would be found. In particular, looking 
at the (7— distances, we see that they match each other 
perfectly. One can deduce that these two models are dif- 
ferent projections in different parameter subspaces of a 
more general EoS parameters vector space (i.e. model). 
One piece of evidence highlighting the equivalence be- 
tween between those scenarios is that the the estimations 
of wq are identical independently of the data used. 

The main difference lies, obviously, in the EoS pivot 
parameter chosen: for CPL it is Wa, whereas Wang takes 
Wo. 5- Although this cannot be falsified, intuition sug- 
gests Wang's pivot is consistently chosen in a "redshift 
sense" , as it is informing us of a dark energy feature 
associated with a redshift value within the observation- 
ally described range. The CPL pivot parameter, on the 
other hand is typically very poorly constrained. Apart 
from the fact that the dark energy evolution in the CPL 
form is largely redshift independent asymptotically at 
high redshifts (lim^^oo {z + 1) = 1), we are question- 
ing ourselves about a physically inaccessible dark energy 
feature, so it is not surprising that for this second dark 
energy parameter we get exceendingly large errors, as 
clearly reflected by our tables. Wang's remedy improves 
the situation substantially; in fact percentual errors on 
Wq and Wo.5 are of the same order of magnitude, whereas 
for Wa they can even be one order of magnitude larger. 



This pattern is common to results coming from real and 
synthetic data. 



2. Dataset choices / combinattons 

In this subsection we look deeper into the results of our 
analysis to gain futher insight on the predictions from dif- 
ferent datasets or their combination. This will allow us to 
draw more definitive conclusions about evolutionary fea- 
tures, to establish how convincing those conclusions are 
for current data, and to forecast what improvements are 
expected from future data. Likewise, these results will 
throw extra light on the effect of different prior choices 
(associated with the still considerable uncertainty in the 
determination of VLm and fif,). 

All real data SNela cases are Icr compatible with 
ACDM. However, BAO data favour phantom models 
quite manifestly as their best fits are concerned, and 
the consequence is that for the SNeIa-|-CMB combina- 
tion ACDM typically is marginally excluded at Icr, i.e. 
the case = — 1 lies outside the Icr boundary but very 
close to it. Arguably, the phantomizing effect of BAO 
data is compensated by their large errors, and possible 
narrower error bands would yield a more decisive exlusion 
of ACDM. 

Mock data, in contrast with real data, yield a picture 
in which ACDM is marginally excluded at Icr by SNela 
data, whereas in the BAO and SNela-l-BAO cases the ex- 
clusion rate is considerable both because our mock BAO 
data favour more phantom models, and because the error 
bands are narrower. 

On the other hand, dark energy subdominance at early 
times, i.e. Woo = wo + u^a < is in general not guaran- 
teed, and there are many cases (both for real and mock 
data) for which that situation is allowed at Icr. The main 
reason for this is that Wa is rather poorly constrained. 
Specifically, for the SNeIa-|-BAO combination from real 
data and at the level of confidence we have just men- 
tioned one can never guarantee Wao < 0, whereas this 
is possible for some mock data cases. Another trend we 
observe is that the lowest values of Woo correspond to 
^l„i — 0.260, the largest value considered. 

A related result is that for SNela the larger the 
lower wq] whereas for BAO the contrary happens and 
drives to the same behaviour for SNela+BAO (of course, 
unless we specify it otherwise, the same pattern is ob- 
served for real and mock data) . This is as well a reflection 
of the tension between the two datasets. In constrast, in- 
creasing f2„i induces always a decrease in Wa- As of fif, 
we can say that when it decreases, so does wq, but the 
opposite happens to Wa- 

In general, if and Sl^ are flxed with independent 
observations, indeterminacies in their values can affect 
our conclusions about the current and early values of the 
dark energy parameter wiyz). The sensitivity to f2m is 
quite noticeable, and becomes particularly manifest for 
the Wa parameter in the CPL case, whereas the situation 
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is not so dramatic in Wang's case. Note that even though 
Wang's parametrization did not try to heal correlation 
between il^ and dark energy parameters, it comes out 
as a valuable consequence. This parametrization throws 
another interesting result: for real data and cither fi,„ < 
0.260 or rib < 0.044 it can be ascertained at la that 
acceleration has increased recently (w has become more 
negative), whereas that situation occurs for mock data 
in all cases. 



C. Main conclusions 

A number of relevant issues about dark energy have 
been explored using SNela and BAO observational data 
within the framework of the most popular representation 
of its possible evolution: the so called CPL model. The 
study is complemented by an in parallel consideration of 
a Wang's reformulation of the former, as this is an ad- 
vantageous alternative which minorates the correlation of 
parameters. Our main goals were to study whether there 
is tension between these two datasets, which are so far 
two of the most worthy tools to explore dark energy, and 
which are anticipated to play an even more preminent 
role in the future, particularly due to the spectacular ad- 
vances expected in BAO astronomical data. 

We have confirmed that tension is present in current 
data, and we have shown that, quite likely, it will be 
present in future data. Regarding this aspect of our work, 
it must be stressed that ours are the first time simula- 
tions of BAO data from both the radial and transverse 
directions in the literature. 

In addition, we have paid considerable attention to the 
influence of indeterminacies in the values of flm and Qi, 



on the tension issue and other cosmodynamical questions 
such as whether dark energy is phantom-like, whether 
there has been an speed up in the acceleration of the uni- 
verse, or whether dark energy could have perhaps been 
the dominant component at early times. Our main con- 
clusions in this respect is that for most priors, and for 
the SNela-fBAO combination best fits, the Universe is 
currently phantom-like, its dark energy EOS parameter 
has become more negative recently, and was dark matter 
dominated at early times. 

Complementary conclusions are that future BAO data 
will improve constraints considerably making them far 
tighter, and Wang's parametrization is a indeed a very 
good way to reformulate the CPL model and could be- 
come the dark energy parametrization of preferred use. 
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Table VI. CPL model and present data. Column. 1; prior choice; columns 2-3-4: best fit values of wq and Wa for different 
datasets combinations; columns. 5-6-7-8: calculated as described in Sec. \V\ with respect to different datasets model. 



Priors 






Best fit parameters {wo, Wa) 






(7 — 


distances 




SN 


BAO 


olN +B AU 


tSN 

d„ 


,BAO 
da 


da 


^BAO-SN 


(0.214, 0.0444) 


-0.965 ±0.128, 


0.540 ± 0.736 


— 1.476 ± 0.713, —0.645 ± 6.955 


— 1.129 ± 0.119, 1.278 ±0.697 


1.04 


> 3 


i.yo 


^ o 


(0.2o7, 0.0444) 


-0.983 ±0.137, 


0.367 ±0.804 


1 I'yl 1 A 'TOO O COT 1 '7 OCC 

— 1.171 ± 0.738, —2.687 ± 7.255 


1 in^y 1 A100 Ani'7 i nTcc 

— 1.107 ± 0.128, 0.917 ± 0.765 


0.62 


> 3 


0.88 


> 3 


(0.260,0.0444) 


-0.999 ±0.148, 


0.140 ±0.885 


-0.880 ± 0.763, -4.666 ± 7.568 


-1.084 ± 0.138, 0.503 ± 0.848 


0.30 


> 3 


0.21 


> 3 


(0.214,0.0405) 


-0.965 ±0.128, 


0.540 ± 0.736 


-1.629 ± 0.725, 0.046 ± 7.099 


-1.153 ±0.119, 1.389 ±0.696 


1.27 


> 3 


1.98 


> 3 


(0.237,0.0405) 


-0.983 ±0.137, 


0.367 ±0.804 


-1.326 ± 0.752, -2.009 ± 7.423 


-1.132 ±0.128, 1.030 ±0.764 


0.81 


> 3 


0.94 


> 3 


(0.260,0.0405) 


-0.999 ±0.148, 


0.140 ±0.885 


-1.036 ± 0.779, -4.006 ± 7.762 


-1.109 ± 0.138, 0.620 ± 0.846 


0.46 


> 3 


0.34 


> 3 


(0.214,0.0366) 


-0.965 ±0.128, 


0.540 ± 0.736 


-1.799 ± 0.740, 0.805 ± 7.272 


-1.180 ±0.119, 1.514 ±0.694 


1.53 


> 3 


2.03 


> 3 


(0.237,0.0366) 


-0.983 ±0.137, 


0.367 ±0.804 


-1.498 ± 0.768, -1.269 ± 7.624 


-1.159 ± 0.128, 1.158 ±0.762 


1.05 


> 3 


1.03 


> 3 


(0.260,0.0366) 


-0.999 ±0.148, 


0.140 ±0.885 


-1.210 ± 0.798, -3.289 ± 7.997 


-1.137 ± 0.138, 0.752 ± 0.844 


0.66 


> 3 


0.52 


> 3 



Table VII. Wang model and present data. Column. 1: prior choice; columns 2-3-4: best fit values of wq and wq.s for different 
datasets combinations; columns. 5-6-7-8: calculated as described in Sec. (|v| with respect to different datasets model. 



Priors 






Best fit parameters (wo, wq.s) 








a - 


- distances 




{n,„,ni,) 


SN 




BAO 


SN±BAO 






^ACDM 


^BAO-SN 


(0.214,0.0444) 


-0.965 ±0.128, - 


-0.785 ±0.132 


-1.476 ±0.713, 


-1.691 ± 1.640 


-1.129 ±0.119, 


-0.703 ±0.128 


1.04 


> 3 


1.96 


> 3 


(0.237,0.0444) 


-0.983 ±0.137, - 


-0.861 ±0.146 


-1.171 ±0.738, 


-2.067 ± 1.717 


-1.107 ±0.128, 


-0.802 ±0.143 


0.62 


> 3 


0.88 


> 3 


(0.260,0.0444) 


-0.999 ±0.148, - 


-0.953 ±0.164 


-0.880 ±0.763, 


-2.436 ± 1.797 


-1.084 ±0.138, 


-0.916 ±0.161 


0.30 


> 3 


0.21 


> 3 


(0.214,0.0405) 


-0.965 ±0.128, - 


-0.785 ±0.132 


-1.629 ±0.725, 


-1.613 ± 1.677 


-1.153 ±0.119, 


-0.690 ±0.127 


1.27 


> 3 


1.98 


> 3 


(0.237,0.0405) 


-0.983 ±0.137, - 


-0.861 ±0.146 


-1.326 ±0.752, 


-1.996 ± 1.759 


-1.132 ±0.128, 


-0.788 ±0.142 


0.81 


> 3 


0.94 


> 3 


(0.260,0.0405) 


-0.999 ±0.148, - 


-0.953 ±0.164 


-1.036 ±0.779, 


-2.372 ± 1.847 


-1.109 ±0.138, 


-0.902 ± 0.160 


0.46 


> 3 


0.34 


> 3 


(0.214,0.0366) 


-0.965 ±0.128, - 


-0.785 ±0.132 


-1.799 ±0.740, 


-1.530 ± 1.720 


-1.180 ±0.119, 


-0.675 ±0.127 


1.53 


> 3 


2.03 


> 3 


(0.237,0.0366) 


-0.983 ±0.137, - 


-0.861 ±0.146 


-1.498 ±0.768, 


-1.921 ± 1.811 


-1.159 ±0.128, 


-0.773 ±0.142 


1.05 


> 3 


1.03 


> 3 


(0.260,0.0366) 


-0.999 ±0.148, - 


-0.953 ±0.164 


-1.210 ±0.798, 


-2.306 ± 1.907 


-1.137 ±0.138, 


-0.886 ±0.160 


> 3 


0.52 


> 3 





20 



Table VIII. CPL model and mock data. Column. 1: prior choice; columns 2-3-4: best fit values of wo and Wa for different 
dataset combinations; columns. 5-6-7-8: da calculated as described in Sec. (|V|. 



CMB - oriented 



Priors 




Best fit parameters (wo,Wa) 






tj — 


distances 




(f2„,Sli,) 


SN 


BAO 


SN±BAO 


dT 




^ACDM 


^BAO-SN 


(0.214,0.0444) 


-1.058 ± 0.057, 0.565 ± 0.274 


-1.383 ± 0.015, 1.415 ± 0.055 


-1.347 ± 0.014, 1.314 ± 0.053 


> 3 


2.62 


> 3 


> 3 


(0.237, 0.0444) 


-1.079 ± 0.061, 0.366 ± 0.299 


-1.276 ± 0.016, 0.990 ± 0.064 


-1.256 ± 0.016, 0.932 ± 0.063 


> 3 


0.93 


> 3 


> 3 


(0.260,0.0444) 


-1.096 ± 0.067, 0.100 ± 0.328 


-1.1.53 ± 0.019, 0.489 ± 0.078 


-1.152 ±0.018, 0.476 ±0.076 


> 3 


0.06 


> 3 


> 3 


(0.214,0.0405) 


-1.058 ± 0.057, 0.565 ± 0.274 


-1.422 ± 0.015, 1.460 ± 0.057 


-1.380 ±0.014, 1.343 ±0.055 


> 3 


> 3 


> 3 


> 3 


(0.237,0.0405) 


-1.079 ± 0.061, 0.366 ± 0.299 


-1.313 ± 0.017, 1.015 ± 0.068 


-1.287 ± 0.016, 0.941 ± 0.066 


> 3 


1.44 


> 3 


> 3 


(0.260,0.0405) 


-1.096 ± 0.067, 0.100 ± 0.328 


-1.187 ±0.019, 0.484 ±0.083 


-1.179 ± 0.019, 0.456 ± 0.080 


> 3 


0.11 


> 3 


> 3 


(0.214,0.0366) 


-1.058 ± 0.057, 0.565 ±0.274 


-1.462 ± 0.015, 1.502 ± 0.059 


-1.412 ± 0.015, 1.369 ± 0.0-58 


> 3 


> 3 


> 3 


> 3 


(0.237,0.0366) 


-1.079 ± 0.061, 0.366 ± 0.299 


-1.349 ±0.017, 1.034 ±0.071 


-1.317 ± 0.017, 0.944 ± 0.069 


> 3 


1.97 


> 3 


> 3 


(0.260,0.0366) 


-1.096 ± 0.067, 0.100 ± 0.328 


-1.219 ± 0.020, 0.470 ± 0.089 


-1.205 ± 0.019, 0.427 ± 0.086 


> 3 


0.42 


> 3 


> 3 


BAO — oriented 


Priors 




Best fit parameters (wo,Wa) 






cr — 


distances 




(0„,06) 


SN 


BAO 


SN±BAO 


df 




^ACDM 


^BAO-SN 


(0.214,0.0444) 


-1.019 ± 0.057, 0.368 ± 0.278 


-1.387 ±0.015, 1.324 ±0.058 


-1.344 ± 0.015, 1.206 ± 0.057 


> 3 


> 3 


> 3 


> 3 


(0.237, 0.0444) 


-1.031 ± 0.062, 0.123 ± 0.304 


-1.270 ±0.017, 0.845 ±0.070 


-1.245 ± 0.017, 0.768 ± 0.069 


> 3 


1.36 


> 3 


> 3 


(0.260,0.0444) 


-1.037 ± 0.067, -0.198 ± 0.333 


-1.134 ± 0.020, 0.265 ± 0.088 


-1.127 ± 0.019, 0.233 ± 0.085 


> 3 


0.10 


> 3 


> 3 


(0.214,0.0405) 


-1.019 ± 0.057, 0.368 ±0.278 


-1.425 ± 0.016, 1.360 ± 0.061 


-1.376 ± 0.015, 1.225 ±0.060 


> 3 


> 3 


> 3 


> 3 


(0.237, 0.0405) 


-1.031 ± 0.062, 0.123 ± 0.304 


-1.305 ± 0.018, 0.854 ± 0.074 


-1.273 ± 0.017, 0.761 ± 0.072 


> 3 


1.89 


> 3 


> 3 


(0.260, 0.0405) 


-1.037 ± 0.067, -0.198 ± 0.333 


-1.163 ± 0.021, 0.233 ± 0.094 


-1.149 ± 0.020, 0.186 ± 0.091 


> 3 


0.36 


> 3 


> 3 


(0.214,0.0366) 


-1.019 ± 0.057, 0.368 ± 0.278 


-1.463 ± 0.016, 1.391 ± 0.063 


-1.407 ± 0.016, 1.239 ± 0.062 


> 3 


> 3 


> 3 


> 3 


(0.237,0.0366) 


-1.031 ± 0.062, 0.123 ± 0.304 


-1.339 ±0.019, 0.855 ±0.079 


-1.300 ± 0.018, 0.746 ± 0.077 


> 3 


2.43 


> 3 


> 3 


(0.260,0.0366) 


-1.037 ± 0.067, -0.198 ± 0.333 


-1.190 ± 0.022, 0.186 ± 0.102 


-1.171 ± 0.021, 0.128 ± 0.098 


> 3 


0.79 


> 3 


> 3 


SNela — oriented 


Priors 




Best fit parameters (wo,Wa) 






cr — 


distances 




(fi^.n^) 


SN 


BAO 


SN±BAO 


df 




^ACDM 


^BAO-SN 


(0.214,0.0444) 


-0.850 ± 0.055, 0.064 ± 0.266 


-1.308 ± 0.014, 1.280 ± 0.052 


-1.257 ± 0.013, 1.137 ±0.051 


> 3 


> 3 


> 3 


> 3 


(0.237, 0.0444) 


-0.860 ± 0.059, -0.136 ± 0.289 


-1.200 ± 0.016, 0.867 ± 0.061 


-1.164 ± 0.015, 0.762 ± 0.060 


> 3 


2.39 


> 3 


> 3 


(0.260,0.0444) 


-0.865 ± 0.064, -0.390 ± 0.314 


-1.080 ±0.018, 0.391 ±0.072 


-1.060 ± 0.017, 0.325 ±0.071 


> 3 


0.81 


> 3 


> 3 


(0.214,0.0405) 


-0.850 ± 0.055, 0.064 ± 0.266 


-1.347 ±0.014, 1.326 ±0.054 


-1.289 ± 0.014, 1.166 ± 0.053 


> 3 


> 3 


> 3 


> 3 


(0.237,0.0405) 


-0.860 ± 0.059, -0.136 ± 0.289 


-1.237 ± 0.016, 0.897 ± 0.064 


-1.194 ± 0.015, 0.774 ± 0.062 


> 3 


2.93 


> 3 


> 3 


(0.260,0.0405) 


-0.865 ± 0.064, -0.390 ± 0.314 


-1.113 ±0.018, 0.395 ±0.077 


-1.086 ±0.018, 0.312 ±0.075 


> 3 


1.30 


> 3 


> 3 


(0.214,0.0366) 


-0.850 ± 0.055, 0.064 ± 0.266 


-1.385 ±0.015, 1.370 ±0.056 


-1.320 ± 0.014, 1.192 ±0.055 


> 3 


> 3 


> 3 


> 3 


(0.237,0.0366) 


-0.860 ± 0.059, -0.136 ± 0.289 


-1.273 ± 0.017, 0.921 ± 0.067 


-1.224 ± 0.016, 0.781 ± 0.065 


> 3 


> 3 


> 3 


> 3 


(0.260,0.0366) 


-0.865 ± 0.064, -0.390 ± 0.314 


-1.145 ± 0.019, 0.390 ± 0.081 


-1.112 ± 0.018, 0.291 ± 0.079 


> 3 


1.82 


> 3 


> 3 
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Table IX. Wang model and mock data. Column. 1: prior choice; columns 2-3-4: best fit values of wq and wo. 5 for different 
dataset combinations; columns. 5-6-7-8: da calculated as described in Sec. (|V|. 



CMB - oriented 



Priors 






Best fit parameters (wq, W0.5) 








a - 


- distances 






SN 




BAO 


SN+BAO 


(i™ 




^ACDM 


^BAO-SN 


(0.214,0.0444) 


-1.058 ±0.057, - 


-0.870 ± 0.040 


-1.383 ±0.015, 


-0.911 ± 0.007 


-1.347 ±0.014, 


-0.908 ± 0.007 


> 3 


2.62 


> 3 


> 3 


(0.237, 0.0444) 


-1.079 ±0.061, - 


-0.957 ± 0.045 


-1.276 ± 0.016, 


-0.946 ± 0.008 


-1.256 ±0.016, 


-0.946 ± 0.008 


> 3 


0.93 


> 3 


> 3 


(0.260,0.0444) 


-1.096 ±0.067, - 


1.063 ± 0.050 


-1.153 ±0.019, 


-0.990 ±0.010 


-1.152 ±0.018, 


-0.993 ±0.010 


> 3 


0.06 


> 3 


> 3 


(0.214,0.0405) 


-1.058 ±0.057, - 


-0.870 ± 0.040 


-1.422 ±0.015, 


-0.936 ± 0.007 


-1.380 ±0.014, 


-0.932 ±0.007 


> 3 


> 3 


> 3 


> 3 


(0.237,0.0405) 


-1.079 ±0.061, - 


-0.957 ± 0.045 


-1.313 ±0.017, 


-0.974 ± 0.009 


-1.287 ±0.016, 


-0.973 ± 0.009 


> 3 


1.44 


> 3 


> 3 


(0.260,0.0405) 


-1.096 ±0.067, - 


1.063 ± 0.050 


-1.187 ± 0.019, 


-1.025 ±0.011 


-1.179 ±0.019, 


-1.027 ±0.011 


> 3 


0.11 


> 3 


> 3 


(0.214,0.0366) 


-1.058 ±0.057, - 


-0.870 ± 0.040 


-1.462 ±0.015, 


-0.961 ± 0.007 


-1.412 ±0.015, 


-0.956 ±0.007 


> 3 


> 3 


> 3 


> 3 


(0.237,0.0366) 


-1.079 ±0.061, - 


-0.957 ± 0.045 


-1.349 ±0.017, 


-1.005 ±0.009 


-1.317 ±0.017, 


-1.002 ±0.009 


> 3 


1.97 


> 3 


> 3 


(0.260,0.0366) 


-1.096 ±0.067, - 


1.063 ± 0.050 


-1.219 ± 0.020, 


-1.062 ±0.012 


-1.205 ±0.019, 


-1.062 ±0.012 


> 3 


0.42 


> 3 


> 3 



BAO — oriented 



Priors 






Best fit parameters (wo,wo.5) 








a - 


- distances 




(a„,,06) 


SN 




BAO 


SN+BAO 




dT° 


^ACDM 


^BAO-SN 


(0.214,0.0444) 


-1.019 ±0.057, - 


-0.896 ± 0.041 


-1.387 ± 0.015, 


-0.945 ±0.007 


-1.344 ±0.015, 


-0.942 ±0.007 


> 3 


> 3 


> 3 


> 3 


(0.237, 0.0444) 


-1.031 ±0.062, - 


-0.990 ± 0.046 


-1.270 ±0.017, 


-0.989 ± 0.009 


-1.245 ±0.017, 


-0.988 ± 0.009 


> 3 


1.36 


> 3 


> 3 


(0.260, 0.0444) 


-1.037 ±0.067, - 


1.104 ± 0.051 


-1.134 ± 0.020, 


-1.046 ±0.012 


-1.127 ±0.019, 


-1.049 ±0.012 


> 3 


0.10 


> 3 


> 3 


(0.214,0.0405) 


-1.019 ±0.057, - 


-0.896 ± 0.041 


-1.425 ±0.016, 


-0.972 ±0.008 


-1.376 ±0.015, 


-0.968 ±0.007 


> 3 


> 3 


> 3 


> 3 


(0.237,0.0405) 


-1.031 ±0.062, - 


-0.990 ± 0.046 


-1.305 ±0.018, 


-1.020 ±0.010 


-1.273 ±0.017, 


-1.019 ±0.010 


> 3 


1.89 


> 3 


> 3 


(0.260,0.0405) 


-1.037 ±0.067, - 


0.198 ± 0.333 


-1.163 ±0.021, 


-1.086 ±0.013 


-1.149 ±0.020, 


-1.087 ±0.013 


> 3 


0.36 


> 3 


> 3 


(0.214,0.0366) 


-1.019 ±0.057, - 


-0.896 ± 0.041 


-1.463 ±0.016, 


-0.999 ± 0.063 


-1.407 ±0.016, 


-0.994 ±0.008 


> 3 


> 3 


> 3 


> 3 


(0.237,0.0366) 


-1.031 ±0.062, - 


-0.990 ± 0.046 


-1.339 ±0.019, 


-1.054 ± 0.010 


-1.300 ±0.018, 


-1.051 ±0.010 


> 3 


2.43 


> 3 


> 3 


(0.260,0.0366) 


-1.037 ±0.067, - 


-1.104 ± 0.051 


-1.190 ± 0.022, 


-1.128 ±0.014 


-1.171 ±0.021, 


-1.128 ±0.014 


> 3 


0.79 


> 3 


> 3 



SNela — oriented 



Priors 






Best fit parameters (wo.wo.o) 








a - 


- distances 




(a™,0(,) 


SN 




BAO 


SN+BAO 


df 


dT° 


^ACDM 


^BAO-SN 


(0.214,0.0444) 


-0.8.50 ±0.055, - 


-0.829 ± 0.039 


-1.308 ±0.014, 


-0.882 ±0.006 


-1.257 ±0.013, 


-0.878 ± 0.006 


> 3 


> 3 


> 3 


> 3 


(0.237, 0.0444) 


-0.860 ±0.059, - 


-0.905 ± 0.043 


-1.200 ±0.016, 


-0.911 ± 0.008 


-1.164 ±0.015, 


-0.910 ±0.007 


> 3 


2.39 


> 3 


> 3 


(0.260, 0.0444) 


-0.865 ±0.064, - 


-0.995 ± 0.048 


-1.080 ±0.018, 


-0.949 ± 0.009 


-1.060 ±0.017, 


-0.951 ±0.009 


> 3 


0.81 


> 3 


> 3 


(0.214,0.0405) 


-0.850 ±0.055, - 


-0.829 ± 0.039 


-1.347 ± 0.014, 


-0.905 ± 0.006 


-1.289 ±0.014, 


-0.900 ± 0.006 


> 3 


> 3 


> 3 


> 3 


(0.237,0.0405) 


-0.860 ±0.059, - 


-0.905 ± 0.043 


-1.237 ± 0.016, 


-0.938 ± 0.008 


-1.194 ±0.015, 


-0.936 ±0.008 


> 3 


2.93 


> 3 


> 3 


(0.260,0.0405) 


-0.865 ± 0.064, - 


-0.995 ± 0.048 


-1.113 ±0.018, 


-0.981 ± 0.010 


-1.086 ±0.018, 


-0.982 ±0.010 


> 3 


1.30 


> 3 


> 3 


(0.214,0.0366) 


-0.850 ±0.055, - 


-0.829 ± 0.039 


-1.385 ±0.015, 


-0.929 ±0.007 


-1.320 ±0.014, 


-0.923 ±0.007 


> 3 


> 3 


> 3 


> 3 


(0.237,0.0366) 


-0.860 ±0.059, - 


0.905 ± 0.043 


-1.273 ±0.017, 


-0.966 ± 0.009 


-1.224 ±0.016, 


-0.963 ± 0.008 


> 3 


> 3 


> 3 


> 3 


(0.260,0.0366) 


-0.865 ± 0.064, - 


-0.995 ± 0.048 


-1.145 ±0.019, 


-1.015 ±0.011 


-1.112 ±0.018, 


-1.015 ±0.011 


> 3 


1.82 


> 3 


> 3 



